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Abstract 

We construct a simple analytical model to study the effects of 
cosmic strings on the microwave background radiation. Our model 
is based on counting random multiple impulses inflicted on photon 
trajectories by the string network between the time of recombina- 
tion and today. We construct the temperature auto-correlation func- 
tion and use it to obtain the effective power spectrum index n, the 
rms-quadrupole-normalized amplitude Qrms-ps and the rms temper- 
ature variation smoothed on small angular scales. For the values of 
the scaling solution parameters obtained in Rcfs.llC|],[p| we obtain 
n = 1.14 ± 0.5, Q rms -ps = (4.5 ± 1.5)Gfi and (^-) rms = 5.5G^. 
Demanding consistency of these results with the COBE data leads to 
G/i = (1.7 ± 0.7) x 10 -6 (where /x is the string mass per unit length), 
in good agreement with direct normalizations of jx from observations. 



1 Introduction 

The recent detection of anisotropy^] in the Microwave Background Radia- 
tion (MBR) by the COBE (COsmic Background Explorer) collaboration has 
provided a new powerful experimental probe for testing theoretical cosmo- 
logical models. The DMR (Differential Microwave Radiometer) instrument 
of COBE has provided temperature sky maps leading to the rms sky varia- 
tion yj< (^) 2 > (where AT = T(6>i)-T(6> 2 ), and 6>i — 6> 2 = 60° is the beam 
separation in the COBE experiment) and the rms quadrupole amplitude. A 
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fit of the data to spherical harmonic expansions has also provided the angu- 
lar temperature auto-correlation function C(A9) =< > where 
<> denotes averaging over all directions in the sky, 5T{6) = T{6)— < T > 
and AO = 6 — 9'. This result was then used to obtain the rms-quadrupole- 
normalized amplitude Q r ms-PS an d the index n of the power law fluctuation 
spectrum assumed to be of the form P{k) ~ k n . The published results are: 

AT 

(^r)rms = (1.1 ±0.2) x 1(T 5 (1) 

Qrms-PS = (5.96 ± 0.75) x 10- 6 (2) 

n = 1.1 ±0.5 (3) 

These results have imposed severe constraints on several cosmological models 
for large scale structure formation. Even the standard CDM model with bias 
1.5 < &8 < 2.5 is inconsistent with the COBE results for Hq > 50km/(sec ■ 
Mpc) and is barely consistent for Hq ~ 50km/ (sec ■ Mpc) || fl5[| (Tensor 
mode perturbations have recently been shown however, to make standard 
CDM with specific inflationary models consistent with COBE for a wider 
cosmological parameter region |23[| ), 

It is therefore interesting to investigate the consistency of alternative 
models with respect to the COBE measurments. The natural alternative to 
models based on adiabatic Gaussian perturbations generated during infla- 
tion are models where the primordial perturbations are created by topolog- 
ical defects like cosmic strings global monopoles or textures. 

Cosmological models based on topological defects have a single free pa- 
rameter v, the scale of symmetry breaking leading to defect formation. Con- 
sistency with large scale structure observations, constrains this parameter 



to Gv 2 ~ 10 _6 |T^][18] where G is Newton's constant. In what follows we 
will concentrate on the case of cosmic strings. 

Previous analytical studies of MBR anisotropics due to cosmic strings [[i"4|l 
were based on the old picture of the cosmic string network evolution [|l6|l and 
therefore focused on the effects of cosmic string loops. Loops however were 
later shown by more detailed simulations || [1C] to be unimportant compared 
to long strings. 

More recently, numerical simulations have been used to investigate the 
MBR predictions of cosmic string models || and comparison of these pre- 
dictions has been made with the COBE datajTj. It was found that for 
Gv 2 ~ 10~ 6 , (v 2 = [i where [i is the mass per unit length of the string) 
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cosmic strings are consistent with the COBE data for a wider range of cos- 
mological parameters than the standard CDM model. The numerical anal- 
ysis that has led to this result however, is rather complicated and currently 
there is still some controversy among the different groups HflO|][]|] about the 
details of the simulations involved in the analysis. In addition, studies based 
on string simulations have necessarily fixed scaling solution parameters and 
therefore, the dependence of the results on these parameters can not be re- 
vealed. These arguments make the construction of an analytical model for 
the study of the effects of topological defects on the MBR, a particularly 
interesting prospect. It is such an analytical model that we are constructing 
in this paper. 

In particular, we use a multiple impulse|l^] approximation to obtain the 
temperature auto-correlation function C(A0) predicted by the string model. 
From C( A9) we obtain the mass per unit length ju of strings consistent with 
COBE and the effective power spectrum index n predicted by the string 
model. 

Our basic assumptions and approximations are the following: 

1. We approximate the photon path from t rec to to by a discrete set of 
N Hubble time-steps ti such that tj+i = 2ij. For z rec — 1400 we have 
JV ~ log 2 [(1400) 3 / 2 ] ~ 16. 

2. At each Hubble time-step the photon beam is affected only by the long 
strings within a horizon distance from the beam. The effects of further 
strings are cancelled due to the compensating scalar field radiation. 

3. The combined effects of all strings present within a horizon distance of 
the photon beam is a linear superposition of the individual effects. [17] 



4. Each string that affects the photon beam induces a temperature vari- 
ation of the formplpg]: 



with 

(3 = k ■ (v s j s x s) 

where k, s and v s are the unit wave- vector, the unit vector along the 
string and the string velocity vector respectively (see Fig. 1). 

The long strings within each horizon volume have random velocities 
positions and orientations. 
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6. The effects of loops are unimportant compared to the effects of long 
strings §0. 



7. Initial temperature inhomogeneities at the last scattering surface are 
assumed negligible compared to those induced by the string network 
at later times. 

We will also use the result that for any function 9\ (9) that takes random 
values as the independent variable varies we have: 

< f(Oi(0)) >e=< f(9i) >0i 

where <> Q implies averaging with respect to a. We will call this for obvious 
reasons the 'ergodic hypothesis'. 

Some of these assumptions are similar to those made in corresponding 
numerical studies. Others (like assumption 2 which is an attempt to take 
into account compensation) are improvements over those of the numerical 
analyses. We comment on the possibility of further improving some of these 
assumptions in section 4. 



2 The Temperature Correlation Function 

We begin with a description of our model. Consider a photon beam emerging 
from the last scattering surface at t\ = t rec . This beam of fixed temperature 
will initially suffer M 'kicks' from the M long strings within the horizon at 
t\. At the Hubble time ti ~ 2t rec , M strings, uncorrelated with the previous 
ones will be within the horizon giving the photon M further 'kicks' and the 
process will continue until the N ~ 16 Hubble time-step corresponding to 
the present time to ( *o — 2 16 i rec ). Therefore the total temperature shift 
in the direction 8 (in what follows we consider fixed ip and omit it when 
defining direction unless otherwise needed) may be written as: 

<- T N M 

n=l m=l 

where f3 mn (9) corresponds to the mth string at the nth Hubble time-step. 
Taking v™ n and s mn to be unit vectors along the string velocity and string 
length directions we may write f3 mn {9) = v s j s (k(9)-Rf n (9)) where Rf n {8) = 
v mn x s mn is a unit vector that varies randomly with m and n. 
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It is instructive for what follows to obtain $f-(0) averaged over all direc- 
tions. Denning £ = 4ttG/j,v s j s we have: 

< S -f(9) >=aY,<Ho)-R? n (o) >=z^E / d cos ed<p cos e^e,^ 

n,m m,n 

where 9± is the angle between Ft™* 1 (9) and k{9) (see Fig. 3) and since R™ n (9) 
varies randomly with m, n and with 9, ip on angular scales larger than 9 rec , 
we may use the ergodic hypothesis to obtain: 

ST £~ f 

< -jr(9) >= dcos9 1 dip 1 cos6q = 

m,n 

which is also the naively expected result. The same result could have been 
obtained without using the ergodic hypothesis by simply taking the large 
M x N limit and performing the sum over m,n before the 9, ip integration. 

Our goal is to investigate correlations of fluctuations :4^(0)4j?(0'). Using 
(4) we have: 

rrp srrp N M 

T {9) — {9') = (4vr^) 2 £ £ P™(9)p m ' n ' (9') (5) 

n,n' m,m' 

Define now A9 = 9 — 9' and focus on the case when A9 = 9 P where 9 P is the 
angular size of the horizon at the Hubble time-step t p (1 < p < 16). Long 
strings present at time t will inflict 'kicks' on two photon beams separated 
by A9 = 9 p that are uncorrelated for t < t p but are equal for t > t p when 
A9 is within the horizon scale t (see Fig. 2). Therefore 

R{ k (9) = R( k '(9') iff k > p, j = f, k = k' 

R{ k (9) ^R( k '{9') otherwise 

where ^ here means 'not equal and also uncorrelated'. We may therefore 
split the sum (5) in two parts consisting of correlated and uncorrelated 
products respectively. 

rrp crp N M 

— {9) — {9')=e[ E E Cm-K n (0))(k9')-Rf n '(9>))+ 

n=n'=pm=m'=l 

+ E E Cm ■ knmw) ■ £r' n V))i = m + s 2 ] 

n,n' m,m' 
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where Si refers to the terms of the sum that are correlated on a scale A9 
while S 2 refers to the uncorrelated terms. Averaging Si over all directions 
we obtain 

<Si(A0) >=-Ly / d cos 9d<p cos 9™ n (9, (p)cos(9™ n (9, ^) + A9) 
Att — ' / 

^ tM Qmn u 

where 9™ n is the angle between k{9) and R™ n (9) (see Fig. 3). Since by the 
construction of the model 9\ varies randomly from one correlated patch to 
another we may use the ergodic hypothesis and replace the average over all 
directions with an average over 9\. It is then easy to see that 

< S!(A0) >= £^2jV cor (A0) 

where N cor (A9) = M(N — p(A9)) is the number of terms in Si. For S2 we 
have: 

<S 2 (A0) >=-!- I dcos8dipJ2 cos 07™ coa(0f n ' + AO) = 

Therefore we may write 

< ^(*)fV) >= ^N cor (A8) cos(A0) (6) 

— 1/2 

Using the relations t p = 2 p t rec and A9 = 9t v = z p (for Q Q = 1) where z p 
is the redshift at t p , it is straightforward to show that 

Aft 

p(A8) = 3log 2 (— ) (7) 

"rec 

Using (6) and (7) we obtain 

C(A8) ^< 5 ^(9) 5 ^(9') >= ^M(N - 31og 2 (^)) cos(A0) 

Our assumptions for the construction of C(A9) clearly break down for A9 < 
9 rec ~ 1° (since we must have t p > t rec ) and for A9 > ^ (since we can 
not have N cor < 0). However, the region of validity of (8) may be easily 
extrapolated by slightly shifting A9 by 9 rec ~ 1° to A9 + 9 rec in the log 
(thus extrapolating C(A9) to A9 = 0) and keeping C(A9) = for A9 > \ 
as is physically expected since N cor goes to at large angular separations. 
We may now write the extrapolated auto-correlation function as: 
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C(A9) = ^-M(N - 31og 2 (l + ^-)) cos(A6>) < A9 < ^ 



C(A9) = — < A0 < vr 



(8) 



The validity of this result for C(A6) has also been verified^] by using 
the model introduced here to calculate directly the rms sky variation and 
showing that the result is identical to the one obtained below from (8) and 
(14) (see section 3). In what follows we will use (8) to compare COBE's 
data with the string model predictions. A typical value for M, the number 
of long strings per horizon volume, is M = 10 while for z rec = 1400 we have 
N ~ 16. 

3 Predictions of the Model 

Consider an expansion of the temperature pattern on the celestial sphere in 
spherical harmonics 

Im 

Defining C/ =< |a™| 2 > it may be shown using the addition theorem that 
C(A9) ^< 5 -f(9) d -f(9>) >= i- £(2Z + l)qP,(cos(A0) 

Assuming now a power spectrum of perturbations of the form P(k) ~ k n it 



may be shown [ 15 1 that 



Ul ° 2 r(/+5_«)r(^) [) 

where C2 = ^-(Qrms-Ps) 2 , Qrms-PS being the rms-quadrupole-normalized 
amplitude 0. 

Since our model predicts C(A9), C\ may be obtained for any I using 
the orthogonality relations for the Legendre functions. We may then find n 
and Qrms-PS that give the best fit to (9). Using the symbol-manipulating 
package Mathematical^ we calculated C/ for I = 3... 30. We excluded the 
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quadrupole C2 from our fit (as was done with the COBE data) and con- 
sidered harmonics up to l max = 30 to account for the small angle cutoff of 
COBE due to the finite antena beam size. Our results are rather insensitive 
to the specific cutoff l max for l max > 15. 
Minimizing the sum 

30 p/^ _j_ n^lN-p/^N 

h (Q ~ C2 F0TWP5) ) (10) 

with respect to C2 and n we obtain: 

ra = 1.14 ±0.24 (11) 



Qrms-ps = \I^C 2 = (45±1.5)(GM)e(u«7.)o.i5V^io x lO" 6 (12) 

as the parameters giving the best fit, where (G(i)q, (v s 7^)0.15 and M\q are 
the corresponding quantities in units of 10 -6 , 0.15 and 10 respectively. The 
error bars were y/ < (ni — n)' 2 >i and similarly for C2, where n/(or C\) is 
obtained by equating to the Zth term of the sum (10) while fixing C2(or 
n) to its best fit value. 

In Fig. 4 we show a plot of C\ vs I obtained from C(A9) of (8) (con- 
tinous thin line), superimposed with C/ obtained from (9) for the best fit 
values n = n and Q rm s-PS = Qrms-PS (dashed thick line). The corre- 
sponding results coming from the COBE data are given by (3) and (2) The 
agreement between (11) and (3) is not too surprising since it is well known 
that topological defects with a scaling solution naturally lead to a slightly 
tilted scale invariant Harrisson-Zeldovich spectrum on large angular scales. 
Clearly however, the confirmed prediction of (11) provides a good test show- 
ing that our simple analytic model has fairly robust and realistic features. 

Comparison of the observational result (2) with the prediction (12) leads 
to an estimate of G/u. According to the recent simulations of Ref.[[|, prefered 
values of the scaling solution parameters on the horizon scale are Vg'js ~ 0.15 
and M ~ 10. Using these values we obtain 

Gfi = (1.3 ± 0.5) x 10~ 6 (13) 

This result is in good agreement with the result obtained numerically in 
Ref.§. 

Another way to obtain an estimate for Gfj, is to compare the measured 

{^-)rms = \J < (-^ ^- > T ; r ^ 2 - > ) 2 ~ with the corresponding value predicted by 
our model. 
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The rms temperature variation may be expressed in terms of the temper- 
ature auto-correlation function for a two beam experiment by the relation: 

AT 

{-jrtms = 2(C(0) - C(a)) (14) 

where a is the beam separation angle (60° for COBE). Since however (^f-)1 ms 
for COBE was obtained with a 3.2° Gaussian beam smoothing we must in- 
troduce a similar smoothing before we compare with the COBE result. Thus, 
we first expand C(0) in spherical harmonics and then reconstruct it using 

the window function W(l) = e 2 ( 17 - 8 ) 2 as was done with the COBE data for 
the construction of C(A9). This leads to 

-, 30 

C(A0) = — ]T(2Z + l)QP z (cos(A0))Vy(/) 2 

1=3 

We now use the reconstructed C(9) in (14) to obtain 
AT / 

(y)™ s = V2.4£ 2 M/3 = 0.55(GM)6(tvy,)o.i5V^u x 1(T 5 

As we did for the rms quadrupole variation we may compare the predicted 
rms variation with the observational result (1) of COBE to obtain the value 
of Gfj, that makes our model consistent with the COBE data. For (v s j s ) = 
0.15 and M = 10 the resulting value for Gfj, is: 

Gfi = (2.0 ± 0.5) x 10~ 6 (15) 

consistent with our corresponding result from Q rm s-PS- An estimate of Gfj 
may be obtained by averaging (13) and (15) and extending the error bars 
to cover both (13) and (15) 

Gfj, = (1.7 ± 0.7) x 10~ 6 (16) 

The interesting consequenses of this result will be discused in the next sec- 
tion. 

4 Discussion 

The main result of this letter is that the COBE data constrain the sin- 
gle free parameter Gfj, of the cosmic string model to be in the range given 
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by (16). This range is in good agreement with normalizations of G[i from 
galaxy |lg|l |L9|l and large scale structure formation || [11]. Therefore the cos- 
mic string model for structure formation remains a viable model, consistent 
with the COBE data. 

Recent studies have used large scale structure observations to derive 
the dependence of G[i on the the bias factor b%, defined such that the rms 
fluctuation in a sphere of radius 8/i _1 Mpc is l/6g Q ||] . In Fig. 5 we plot the 
allowed range of b% for four different cosmic string models, demanding that 
the value of G\x in Q^] is consistent with (16). For comparison we also plot 
the allowed range of 6g for the standard CDM model, with cosmological 
parameters chosen to give the most reasonable value for it [15]f|. Our 
estimate for the allowed range of the bias factor in the string model is slightly 
lower than that of Ref.0] but is clearly consistent with it. Since a value of 
b$ in the range 1 to 3 is generally well motivated physically, Fig. 5 suggests 
that both the string model and standard CDM (with h = 0.5) are consistent 
with the COBE data. 

Even though our model is based on fairly simple assumptions it has 
provided results that are not only self-consistent but also in good agreement 
with much more complicated numerical analyses. This fact indicates that 
our assumptions even though simple are fairly realistic and economical. In 
spite of that, there is certainly room for improvement. One interesting 
impovement that could be implemented in a straightforward way is the 
consideration of the effects of loops. Even though we do not expect loops 
to affect the shape of the auto-correlation function on the angular scales 
considered here they may introduce an overall multiplicative factor slightly 
larger than one, thus reducing the predicted value of G/i by the same factor. 
Another improvement that would modify our results in a similar way is the 
consideration of string induced perturbations on the last scattering surface 
by strings present in the final stages of the radiation era. Finally, it would 
be interesting to include the effects of compensation in a more realistic way. 
For example, instead of introducing an abrupt cutoff of the deficit angle on 
the horizon scale we could smoothly reduce it to zero as the horizon scale is 
aproached. The effect of this modification would increase the predicted value 
of Gfi thus tending to cancel the effect of the previous two modifications. 

It is straightforward to apply the analysis presented here to other seed 
based cosmological models including the global monopole and texture mod- 
els. Work in that direction is currently in progress. 
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Figure Captions 

Figure 1 : 

The geometry of the vectors v, s and k. 

Figure 2 : 

Two photon beam paths in directions 9 and 9' from t rec to to- The horizon 
in three Hubble time-steps is also shown. The effects of strings during the 
first two Hubble time-steps are uncorrelated for the particular angular beam 
separation shown. 

Figure 3 : 

The geometry of the vectors k (9) and R r [ in (9) in the term £i(A0). 
Figure 4 : 

The multipole coefficients C\ (continous thin line) superimposed with the 
best fit of (9) (dashed thick line), obtained by using n = 1.14. 

Figure 5 : 

The allowed range of the bias factor bs for four cosmic string based models 
obtained using our result for G/j, and the results of Ref. jjjf|]. The following 
cases are shown: 

a) Strings + CDM, h=l and h=0.5 

b) Strings + HDM, h=l and h=0.5 

The allowed range for standard CDM with h = 0.5 and = 0.1 (the most 
consistent case with COBE's data) is also shown for comparison. 
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